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, Abstract 

We present a quantum mechanical study of the bound states of a neutral scalar particle on 
, the curved background described by the Reissner-Nordstrom (RN) spacetime that corresponds 

to a naked singularity. We show that there occurs both the metastable, i.e. decaying bound 
^ ' . states and the stable ones. The corresponding energy spectra are calculated in the leading 

WKB approximation and the comparison with the RN black hole quasinormal mode spectrum 
$—( ' is made. The metastable bound states on the naked singularity background turn out to be more 

long-living than the quasinormal modes for the black hole with the same mass. 
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1 Introduction 

In general relativity, the probe with test particles and waves is an efficient tool to 
ly-^ I study the properties of spacetime near gravitating mass [THO]. When using quantum 

test particles, one observes some features that are unexpected from the viewpoint of 
classical theory. In particular, it is well established that quantum test particles cannot 
form stable stationary bound states in the field of the Schwarzschild and Reissner- 
^ ! Nordstrom (RN) black holes [7j. Instead, there exist decaying quasistationary states 

. known as the quasinormal modes (see Refs. [8l[9] for review). 

However, if in the RN setup the gravitating mass M and the electric charge Q 
obey the superextremality condition, Q > M, the space-time geometry corresponds 
not to a black hole but to a naked singularity. At the classical level, this results in 
a repulsive potential barrier near the origin that prevents the absorbtion of neutral 
test particles [T014T2] . Therefore one may anticipate that the bound states are stable 
at the quantum level as well. 

The main goal of this article is to verify the above proposition and to derive the 
energy spectrum of the bound states of a quantum test particle on the naked RN 
background. 
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Probing the naked RN singularity, one faces two key problems. The first one 
concerns physical relevance. According to the cosmic censorship conjecture by Penrose 
[T3] . any singularity originated from gravitational collapse must be hidden inside 
an event horizon. Besides, since the RN singularity is electrically charged, it may 
be neutralized by spontaneous pair creation provided that charges are sufficiently 
high [H] . Nevertheless, one cannot exclude scenarios in which the RN geometry with 
Q > M is an exterior solution connected with a regular interior solution with matter, 
for instance, the Friedmann solution in the model of friedmon [ISj. Such models 
can be traced back to Dirac's suggestion to describe leptons as charged shells [IB|. 
Recall that it is Q/M = 2.0 ■ 10^^ for the electron, so that the exterior geometry must 
correspond to the naked singularity. 

The second problem is the non- uniqueness of the time evolution for test particles 
and waves on the naked RN background [Hl[I71[TB]. This means one has to specify an 
additional boundary condition at the singularity to obtain a fully unique dynamics. 
In the present work, we overcome this shortcoming by employing the WKB approach, 
which turns out to be insensible to behavior in the vicinity of the singularity. 

The work uses a simple quantum mechanical model. We consider a neutral scalar 
particle on the background geometry; the word "particle" may refer equally to the 
test scalar field perturbation. Both massless and massive cases are studied. This 
investigation expands our previous work \T^, which was restricted to the stable bound 
states of a massive particle. 

The plan of the article is as follows. In Section [2] we examine how particle wave 
functions behave in the vicinity of the singularity and advocate the WKB approxima- 
tion. In Section [3] the effective potential is analyzed to distinguish between stable and 
decaying bound states. Section H] is devoted to the calculation of the energy spectra 
through the WKB method. The energies obtained are then discussed and compared 
with the known quasinormal mode spectra for the RN black holes in Section [51 Fi- 
nally, Section IH] presents our conclusion. 

2 Behavior in the vicinity of a RN singularity 

To start with, we establish the behavior of a particle wave function in the vicinity 
of a RN singularity and study how it is accommodated in the WKB approximation. 
The external RN geometry is described by the metric (in units with G = h = c = 1) 



For concreteness, we assume that Q > 0. 

The Klein-Gordon equation for a neutral scalar particle of mass m on this back- 
ground 



ds"^ = Fdt^ - F-^dr'^ - r^{de^ + sin^ Odcj)^) , F = 1 



2M g2 



(1) 




(2) 
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has stationary state solutions, which we decompose into the spherical harmonics 

v^(t,r,^,0) =exp(-^)F,„(^,0)^w(r). (3) 

Then the radial wave function satisfies the equation 
cPu 1 



rfr2 ^ F 



1 ( 2 Q^-M^ \ 2 ^(^ + 1) 
' oj — m 



u = (4) 



F \" , ^2 

where Z = 0, 1, 2, ... is the orbital quantum number. In principle, this can be reduced 
to the confluent Heun equation [5D], upon substituting u{r) = exp^y/m"^ — u'^r){M — 
^/M^TQ2_r)»(^M+^/W^^-r)(^y{r) with a and /3 given by bulky expressions 
in terms of Q, M and m. The so-obtained confluent Heun equation has singularities 
at r = M ± a/ — Q'^ (black-hole horizons) and r = oo. However, in the super- 
extreme case under consideration these singularity are shifted away from the real line 
because — < 0. As a consequence, the reduction to the confluent Heun form 
does not simplify the problem. 

It is more useful to note that the radial equation holds invariant under rescaling 

Q aQ, M aM, m — m, u -> — w, R ^ ar (5) 

a a 

and thus one of the three parameters Q, M and m can be fixed arbitrarily. 
Close to the naked singularity at r = 0, the radial equation simplifies 

^^-^±l(i±ipl^. = 0. r^O. (6) 

If we now adopt the leading WKB approximation, choosing the radial wave function 
in the form u{r) oc exp(2S'(r)), then in the super-extreme case we have S"^(r) < 
for small r. Hence, the naked singularity is located in the classically forbidden 
region. As for the bound states we are interested in, this means that the wave 
function asymptotics is dominated by the term which decreases exponentially as r 
approaches 0. Note that for small r the condition of the WKB-method applicability 
\S"{r)/S'\r) \ < 1 is still satisfied. 

Nevertheless, the exact behavior of the wave function at r — )■ remains undeter- 
mined. The reason is that in a pair of linearly independent solutions to Eq. both 
solutions are locally normalized at r — t- and thus none of them can be discarded. 
As seen from Eq. ([6]), such a pair may easily be composed by taking a solution with 
u{0) = and another one with u'{r)\r=o = 0. In the general case, one should consider 
a linear combination of these solutions and employ the mixed boundary condition 
{u'{r) + au{r))\r=o = with an arbitrary real parameter a. Physically, the depen- 
dence of the dynamics on this parameter indicates that the naked singularity carries 
some additional degrees of freedom beyond those contained in the metric. It is said |S] 
that the RN singularity has " hair" . 



3 



Mathematically, the ambiguity in choosing the boundary condition stems from the 
non-uniqueness of the self-adjoint extension to the wave operator on the naked RN 
background [5lll7[ll8j . As a possible resolution, the Friedrichs extension was suggested 
which implies the choice a = and is naturally connected with the quadratic form of 
energy (TTJIIB]- It should be stressed that the application of the WKB method will 
permit us to bypass this problem at all because the method does not uses the value 
of the wave function in the origin. 

In the end of this Section it is worthy to compare the above picture in terms of the 
areal radial coordinate r with that based on the usual tortoise coordinate x defined 
by dr = Fdx. Adjusting the integration constant, one can always put the limit of 
r — )• in correspondence with a; — )■ 0. Then in place of the radial equation (jl]) one 
obtains 



^2 



(7) 



where = u{r)/ y^F{r). In the limit of x — )■ one has x ~ r^/(3Q^) and this 

equation reduces to the Schrodinger-type one with the inverse quadratic potential 

Then the textbook analysis [21] shows that the particle cannot fall on the origin 
because the numerical coefficient of the potential in Eq. (|8]) does not exceed its crit- 
ical value 1/4. This is in agreement with our previous conclusion that the naked 
singularity is located in the classically forbidden region. 

However, Eqs. ([7]) and ([S]) contain the singularity at x = such that the condition 
of the WKB-method applicability breaks down for small x. Therefore, to develop the 
WKB approximation in the subsequent sections, we employ the initial radial equation 
(HI) in terms of the areal radial coordinate r. 



3 Effective potential 

Now we shall establish possible types of the particle bound states on the naked RN 
background. To that end, let us analyze the shape of the effective potential V{r) 
introduced upon rewriting the radial equation (j4]) in shorthand notations 

^ + jrA-'- V(r)) u = 0, (9) 

where 

, 2Mm^ QW + A 2AM Q^A-M^ + Q^ , , 

= ^ ^ + -2 ^ + (10) 
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and we designated A = /(/ + 1). Note that V{r) is positive-defined because F > for 
all positive r. 

For a given particle energy co, the type of the bound state is regulated by the 
pattern of classically allowed and forbidden regions. From Eq. one concludes that 
the classically allowed (forbidden) region is to be defined as a region in which the 
condition u'^ -V{r) > (w^ - V{r) < 0) holds. 

Instead of solving those inequalities explicitly, we are looking for local extrema of 
V{r). If V{r) has a minimum at r = rmin and no other extrema, then for u close 
enough to \fV{j\^^ the classically allowed region consists of a single line segment and 
there stable bound states may exist. In this case we call V{r) the single-well potential. 
If Vir) has two or more extrema, then for certain u the classically allowed region must 
be split into two subregions separated by the classically forbidden region. Because 
of the tunneling between the subregions, the bound states are now metastable, i.e. 
decaying. The corresponding potential will be referred to as the barrier-shaped one. 

We are now in the position to deduce what type of the bound states occurs de- 
pending on the values of the particle mass m and the singularity parameters Q and 
M. Let us first study the simpler case of the massless particle and then turn to the 
massive case. 



3.1 Case of m = 

In this case the effective potential simplifies 

A 2AM , Q^A-M^ + g^ 
^W = ;3-^ + • (11) 

Notice that for the S-waves (/ = ^ A = /(/ + 1) = 0) it is monotonically decreasing 
and thus only scattering states are present in the spectrum. 
For / 7^ the formula 

8 + 9A 



M 8 + 8A 

defines a critical charge-to-mass ratio such that for QjM > {Q/M)cr the effective 
potential V{r) has no extrema. If 1 < Q/M < {Q/M)cr, there must be one local 
minimum and one local maximum, so that the potential is barrier-shaped and there 
emerge the metastable bound states. Noticeably, in the limiting case of the extreme 
black hole, Q/M — )■ 1, the local minimum transforms in an event horizon. 

From Eq. f[T^ . one sees that the value {Q/M)cr increases with the growth of A, 
though not exceeding y/Q/S ~ 1.0607. Thus the metastable bound states may only 
exist in the narrow range 1 < Q/M < {Q/M)cr- The plots of V{r) for m = 0, obtained 
with the values Q/M = 1.01 and 1.1 inside and outside this range, respectively, are 
shown in Fig. 1. 
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(a) (b) 

Figure 1: Effective potential V{r) for the particle with m = on the background with (a) M = 1, 
Q = 1.01, (b) M = 1, Q = 1.1. Different lines refer to the orbital numbers Z = 0, 1, 2. 



3.2 Case of m ^ 

We shall perform the same analysis as in the previous subsection but for the massive 
particle. For brevity, we set M = 1 in this subsection; the formulae with arbitrary 
M can readily be obtained by rescaling 

To study the shape of V{r), we find its extrema that amounts to solving the cubic 
equation 

dV 2im? ,09,, , ^ 

— = ^— (r-^ + ar^ + 6r + c) = (13) 

with the coefficients 

The number of real roots to this cubic equation is controlled by the sign of its deter- 
minant defined by 

where 

0? , oh , , 

, = 2(-) -- + C. (16) 

There must be one, two and three real roots when D > 0, D = and D < 
respectively. Since the signs of the coefficients alternate, all these roots are positive 
and thus all local extrema of Vir) are of physical meaning. 

It is convenient to rewrite D in terms of the deviation from the extreme charge 
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value S 
D 



1 



108m^ 



[(85- 1)A^+ (-78m2(5 + 6m2 + 24mV + 8 5) 



+ (24 mH^ + 15 - 84 + 24 + 63 m^6^ + 54 m^5) 

+ (8 m^6^ + 32 m^6 + 48 m^^^ + 32 m^6^ + 8 + 24 m^6'^ + 156 m^6^ 

+ 132m^5) A + 8m^5 + 24mV + 24m^(5^ + 8mV + 108mV] . (17) 

For the super-extreme charge under consideration, we have 6 > 0. Then it becomes 
evident from f|T7j) that the condition D > always holds true for the S-states (A = 0) 
and, hence, the corresponding effective potential has one local extremum (minimum), 
thereby being single-well. Thus, in contrast to the case of m = 0, the particle with 
m 7^ may form stable bound S-states. 

As an illustration, in Fig. 2 we offer plots of V{r) obtained with m = 0.1, using 
Q/M = 1.01 and 1.1. In these plots the potential well where the bound states emerge 
lies below the rest energy level = rn? depicted by the dotted line. 



0.5 1 

V(r) 

0.4- 



0.3 



0.2 



0.1 




0.5 1 
V{v) 

0.4- 



0.3 



0.2 



0.1 




10 



10 



(a) (b) 
Figure 2: The same as Fig. 1 but for the particle with m = 0.1. 

As for the spectrum with / 7^ 0, we have learnt that in the case of m = it contains 
the metastable bound states provided the charge-to-mass ratio is close enough to its 
extreme value Q/M = 1. We now argue that in the case of m 7^ the situation is 
essentially the same. Indeed, for small 6 the asymptotics of D is given by 

A (8m2-A) (A + m^)^ 



D 



+ 0{6), 



108m8 

so that D < if A > 8m^. In practice, the last condition holds for all / 7^ 0, 
because m <^ M = 1 (the particle mass is supposed to be much less than the central 
object mass). The inequality D < means that dV/dr has three zeros and, as a 
consequence, V{r) is the barrier-shaped potential with two minima and one maximum. 
This subcase may include the metastable bound states because one of the minima 
lies above the level of = (see Fig. 2a). 
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On the other hand, if Q/M is far from being extreme, the metastable bound states 
are excluded. In the case of m = 0, the critical value ( !T2|) that provides the necessary 
and sufficient condition for their exclusion was derived. Now, for m 7^ 0, we are able 
to deduce the sufficient but not necessary condition. To do this, we should establish 
for which values of 6 the condition D > holds, which guarantees that the effective 
potential is single-well. 

Upon demanding that the coefficients in front of powers of A in f|T7|) be positive, 
one obtains several restrictions. The strongest one comes from the coefficient of 
the A^-term which becomes positive-defined when 6 > 7/2. Thus we conclude that 
the condition D > holds and, hence, the effective potential is single-well and the 
metastable bound states are excluded, provided 

It should be stressed that this is not an exact critical value and, in fact, D > may 
hold at lesser values of QjM when / is small enough. For example, if Q/M = 1.1 < 
2.121 the effective potential turns out to be single- well for / = 0, 1, 2. In the case of 
/ = it can be seen from Fig. 2b whereas for / = 1, 2 the bottom of the well is located 
in the region of large r outside the plot. 

Let us summarize our findings. For / = 0, the spectrum contains only the scattering 
states, if the particle is massless, and in addition the stable bound states, if m 7^ 0. 
For I 7^ 0, the metastable bound states emerge if the central object charge-to-mass 
ratio is close to its extreme value QjM = 1. When Q/M exceeds a certain critical 
value, the metastable bound states disappear and one observes the scattering states 
and the stable bound states. 



4 Energy spectrum 

In this Section we derive the energy spectrum of the particle (= test scalar field 
perturbation) in the leading WKB approximation. 

Since the background geometry is naked rather than the black- hole one, our calcu- 
lation will differ essentially from typical computations of the black-hole quasinormal 
modes. The latter imply that the solutions to the wave equation go as plane waves 
close to the horizon. Instead, in our case the solutions have to decrease exponentially 
deep into the classically forbidden region in the vicinity of the origin. 

Actually, the barrier-shaped potentials that lead to the metastable bound states 
in our treatment (see Figs, la and 2a) resemble those of Gamow's theory for a-decay 
(Fig. 3). In turn, the single- well potential we obtained with m 7^ and / = resembles 
the attractive Coulomb potential everywhere except for a narrow region close to the 
origin. Therefore we shall adopt the ordinary flat-space WKB formulae |21] with the 
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Langer modification for tlie centrifugal term, A = /(/ + 1) — )■ (/ + 1/2)^, which is usual 
in the Coulomb-like problems [22||23]. 




Figure 3: Potential in the a-decay theory. 



J^J^l-V{r)dr = 7l{n + l/2) (21) 



4.1 Decaying bound states 

Let us start with the metastable bound states. From the wave equation ([2]), it can 
be seen that the quasiclassical momenta k{r) and K(r) for the classically allowed and 
forbidden regions respectively should be defined by 

k{r) = \^^uj^-V{r), K{r) = \^^V{r) - uj''. (20) 
r r 

For the metastable bound states the energy u = ujr + iuj is complex. In the lowest 
WKB approximation its real part is determined by the Bohr-Sommerfeld rule [21] 

F 

where ri 2 are the turning points at the ends of the potential well (see Fig. 3) and n 
denotes the excitation number. 

Further, we derive a Gamow-type formula for the imaginary part of u. Our deriva- 
tion follows the conventional procedure [21]. First, multiplying the wave equation ^ 
by the complex conjugate wave function, we compose the expression u"u — u"u = 
— {u^ — a}^)F~^|up and then integrate it to obtain 

POO 

{u'u-u'u)l_^^ = -AiuRUi F-Vprfr. (22) 

Jo 

Here we took into account that {u'u — u'u) vanishes at r = because u must satisfy 
the boundary condition u'(0) + au{Q) = with real a in order that the Klein-Gordon 
operator be symmetric. 

Next, we resort to the well-known WKB formulae for the wave function asymp- 
totics in the classically allowed region pT] 



u ~ — ^ exp ( i I k{r)dr — in/A] , r > r^, (23) 
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u ~ — COS ( / k{r)dr — 7r/4 ) , ri < r < r2- (24) 



Using the first of tliese formulae, tlie left-liand side of fl22l) is estimated to be equal 
to i|Cip/2. The main contribution to the right-hand side of fl22|) comes from the 
region ri < r < r2- Approximating the squared cosine in this region by 1/2, we 
rewrite (12^ as 



2 - ^1 -I 2F2A;(r) 

Applying the connection formula [21] C2 = Ci exp{J^^ n{r)dr), we end up with the 
requisite expression for the imaginary part of u 



Ui = — exp I —2 / K{r)dr 



r-3 



r-2 



R 



^2 



dr 



ri 



(25) 



This formula as well as the WKB asymptotics fl23|) and flM|) is valid if the barrier is 
nearly impenetrable. It means that the calculated value of uj, which determines the 
decay rate, has to obey the condition ^ ur. 

4.2 Stable bound states 

Now let us turn to the stable bound states occurring in the case of the massive particle. 
If m 7^ 0, the effective potential (fTOj) has the Coulomb tail {V oc — 1/r), so that a 
Balmer-type approximate formula for the real bound-state energies u = ujr can be 
obtained. 

To do this, we first compare the wave equation dH]) for our effective potential and 
the ordinary Schrodinger equation for the Coulomb potential 

^ + 2m[E + -K^^ = 0. (26) 



dr"^ \ r 2mr2 

Since at large r the factor F of the RN metrics ([T]) tends to 1, it becomes evident 
that the term 2Mm'^/r in Eq. corresponds to 2mZe^ /r in Eq. (12^ . Then we can 
define the characteristic radius for our system, = l/{m?M), corresponding to the 
Bohr radius, l/{mZe^), for Eq. (126|1 . Introducing now the dimensionless coordinate 
p — r /tbi the effective potential ( ITO|) with the Langer modification is rewritten as 

2 2M'^m^ M'^m^lQ'^m^ + il + l/2f] 2il + l/2fM^m^ 

V{p) = \ ^— '-^ ^ '-^ 

p p^ p-^ 

M^m\Q^[l + {l + l/2f]-M^} 

The last two terms in this expression contain higher powers of the particle mass m 
which is much less than the central object mass M. Hence, we may neglect these 
terms provided that M is not much larger than the Planck mass (equal to 1 in our 
units). Thus we obtain the truncated potential 

V (n) 2MW MW[gW + (/ + 1/2)2] 

ytruncKP) = m 1 [Zti) 

p p^ 
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serving as a good approximation to V{p) in the classically allowed region. Using the 
same reasoning, we approximate F = 1 — 2M'^m? / p + M'^Q'^m^ / by 1. 

Within the above approximation, the Bohr-Sommerfeld integral (l2Ti) reduces to 




, 2M2m4 M2m4[Q2^2 + (/ + 1/2)2] 
co^ — m2 H 



P 



P' 



m?M 



7r(n + l/2) (29) 



where turning points pi^2 are zeros of the expression under the square root. This 
integral is easily calculated, by applying the formula 



P2 



V(P-P1)(P2 -P) 



dp 



Pi + P2 



VPiP2 



'PI P 

As a result, we get the explicit expression for the stable bound-state energies 

1 1/2 



(30) 



m 



{n + 1/2 + v/(/ + 1/2)2 + g 



2m2 



(31) 



It has essentially the same structure as the Balmer formula in the Coulomb problem, 
assuming that the values of the quantum numbers are high, n, / ^ 1, as usual in the 
WKB method. The only difference is that the coupling constant is now given by the 
product of masses, but not charges. Note that the quantization procedure that starts 
from the classical particle Hamiltonian has lead us to the same expression 



5 Numerical results 

First we examine the spectrum of the metastable bound states for the massless neutral 
scalar particle in the field of naked RN singularity. It should be noticed that since 
the depth and the width of the corresponding potential well are finite, the number of 
the metastable bound states is limited. The energies of all the existing states with 
/ = 1,2,3 calculated according to the WKB formulae and ([2SD with Q = 1.01, 
M = 1 are presented in Table 1. That the states with larger excitation numbers n 
do not exist was fixed by observing that the Bohr-Sommerfeld integral (12T|) cannot 
be saturated to its value 7r(ra -|- 1/2) even at the maximal allowed bound-state energy 
cj^ equal to the peak of the effective potential. As seen from Table 1, for higher / the 
number of the existing bound states and their half-lives increase as the peak of the 
effective potential grows with I (see Fig. la). 

For comparison, in Table 2 we list the first three quasinormal energies (frequencies) 
of the scalar field on the extreme RN black-hole background with Q = M = 1 that 
were computed in Ref. [25] . 

The striking difference between the results in Tables 1 and 2 can be ascribed 
to the fact that the event horizon is absent in the naked singularity case and thus 
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Table 1: Energies of the metastable bound states on the naked RN background calculated with 
g = 1.01, M = 1 andm = 0. 
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V ^ min 


/T/ 
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0.294636 
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Q 
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0.241152 


0.380747 








/ = 


= 2 












^WKB 
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0.423194 


-i 0.164 X 10" 
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0.364334 
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= 1 


0.537046 


-i 0.227 X 10" 


-3 












/ = 


= 3 
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0.673735 
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0.786160 
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-3 
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= 3 


0.884745 


-i 0.523 X 10" 


-2 







Table 2: Quasinormal energies of the extreme RN black hole with Q — M = \ and m = taken 
from Ref. [25] 
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^WKB 
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0.37764 -z 0.08936 





0.375 
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1 


0.34392 - i 0.27828 






n = 
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0.29661 - i 0.48145 
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^WKB 




^/^max 


n = 





0.62609 -^0.08873 





0.625 


n = 


1 


0.60677-^0.26944 






n = 


2 


0.57254 -z 0.45750 
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the metastable bound states and the quasinormal modes are of different nature. As 
was discussed in the previous Section, the setup for the naked RN bound states 
resembles that of Gamow's theory for a-decay. In this setup the most longhving 
states are concentrated near the bottom of the effective potential well uj^ ~ min[y(r)]. 
This can be readily checked by comparing the calculated energies lo with the values 
= min[\^(r)] which are also shown in Table 1. In turn, the black-hole quasinormal 
modes are, roughly speaking, the scattering resonances travelling both to r = and 
to r = cxD. Thus they have much smaller half-lives and survive for energies close to 
the peak of the effective potential uj^ ~ max[V^(r)] (see Table 2). 

For the naked RN singularity, the quasinormal modes of the latter type were 
calculated in the recent work [26]. These authors postulate the Dirichlet boundary 
condition that rules out one of the linearly independent solutions to the Klein-Gordon 
equation. In the present work we do not specify the boundary condition because it is 
not needed in the WKB approximation. Moreover, for the metastable bound states 
we consider the point r = is located deep into the classically forbidden region and 
thus the boundary condition cannot affect these states substantially. 

Table 3: Energies of the metastable bound states on the naked RN background calculated with 
Q = 1.01, M = 1 and m = 0.1. 
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Next we investigate the case of the massive particle. In Table 3 the energies of all 
the states with / = 1, 2, 3 computed using Q = 1.01, M = 1 and m = 0.1 are hsted. 
Also, in Fig. 4 the imaginary part of u is plotted versus the real part of u for different 
values of the particle mass m. From Tables 1, 3 and Fig. 4, we see that the real part 
of u, i.e. the oscillation frequency, grows with increase of m, while the imaginary 
part of u, representing the decay rate, falls down. Interestingly, the same behavior 
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was observed for the quasinormal modes of the scalar field on the sub-extreme RN 
black-hole background in Ref. [27]. 
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Figure 4: Real and imaginary parts of the metastable bound-state energies for n = 0, I = 1 using 
(a) M = 1, Q = 1.01, (b) Af — 1, Q — 1.02. The squares left to right present the results obtained 
with m = 0, 0.05, 0.1, 0.15, 0.2 respectively. 



Now let analyze the spectrum of the stable bound states. In Table 4 we list the 
energies {uj/m)]YXB of the ground and the first two excited S-states (/ = 0) computed 
according to the WKB formula fl3Tl) and also the energies {uj/m)num obtained by direct 
numerical integration of the Klein-Gordon equation (jlj) with the Dirichlet boundary 
condition u{0) = 0. The calculation was made using Q = 1.1, M = 1 and various 
m. From Table 4 we see that the binding energy ujf,ind = — (w — m) increases as the 
particle mass m increases. However, this binding energy is more than two orders of 
magnitude smaller than the metastable bound-state energy ujr calculated with the 
same values of m and M (see Table 3). It means that transitions between levels of 
the stable states due to an external perturbation would have much lower frequencies 
than the those due to the decay of the metastable bound states. 

Table 4: Energies of the stable bound S-states on the naked RN background with Q = 1.1, M = 1. 



m = 0.05 


m = 0.1 


m = 0.15 


{u}/m)wKB {(^/m)num 


{uj/m)wKB {(^/m)num 




n = 0.998757 0.998733 


0.995105 0.994710 


0.989266 0.987068 


n = 1 0.999688 0.998685 


0.998764 0.998701 


0.997257 0.996906 


n = 2 0.999861 0.999860 


0.999449 0.999429 


0.998771 0.998661 



6 Conclusion 

In this work we have studied the states of neutral scalar particles in the field of 
the naked RN singularity. It has been established that their energy spectrum is 
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qualitatively different from the quasinormal mode spectrum for the RN black hole. 
The conditions for the bound states to be formed have been found and the possible 
types of these states have been examined. 

The first possible type is the metastable, i.e. decaying bound states. They occur 
provided that the charge-to- mass ratio for the central object, Q/M, is close to its 
extreme value Q/M = 1. If the particle is massless, this is expressed by the inequality 
1 < {Q/Mf < [8 + 9/(/ + l)]/[8 + 8/(/ + 1)] where / is the orbital number. 

Using the WKB method, wc have calculated the metastable bound-state spectrum 
which shows that there is no continuous transition in energies between the RN naked 
singularity and the extremal RN black hole cases. The metastable bound states 
on the naked RN background with Q/M > 1 have much larger live-times than the 
quasinormal modes ior Q/M — 1. This is not surprising because for Q/M > 1 there 
exists no event horizon to fall on, so that the metastable states can decay only at 
the expense of the particle escape to spatial infinity upon tunneling through the wide 
potential barrier. 

As the second possible type of particle states, we identify the stable bound states. 
It should be stressed that these states have no analog in the case of the RN black hole 
and can only be formed by massive particles in the field of the central object with the 
sufficiently high charge-to-mass ratio. Their energy spectrum, calculated by means 
of the WKB method, has the structure similar to that of the Coulomb spectrum. 

Energy gaps between the stable bound states, i.e. transition frequencies, arc sev- 
eral orders of magnitude smaller than the frequencies for the metastable bound states 
obtained with the same masses of the particle and the central object. Thus the lat- 
ter states are more favored to be ever found in experiment. Nevertheless, the stable 
bound states may manifest themselves in a different way, by making up a scalar con- 
densate around the central object. Then the question arises as to whether such a 
condensate can provide mass needed to screen the RN singularity. Answering this 
question requires a self-consistent model for interacting scalar and gravitational fields 
which may be constructed in a future work. 
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